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ABSTRACT 

We perform helioseismic holography on realistic solar convection simulations 
and compare the observed travel-time perturbations with the expected travel 
times from the horizontal flows in the simulations computed from forward models 
under the assumption of the Born approximation. We demonstrate reasonable 
agreement between the observed and model travel times which reinforces the 
validity of helioseismic holography in the detection of subsurface horizontal flows. 
An assessment is made of the uncertainty of the measured p-mode travel-times, 
from the root-mean-square of the residuals. From the variation of the signal-to- 
noise ratio with depth, we conclude that the helioseismic detection of individual 
flow structures with spatial scales of supergranulation or smaller is not possible 
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for depths below about 5 Mm below the surface over time scales less than a day. 
The travel-time noise estimated from these simulations appears to be similar to 
noise in measurements made using solar observations. We therefore suggest that 
similar limitations exist regarding the detection of analogous subsurface flows in 
the Sun. A study of the depth dependence of the contribution to the travel- 
time perturbations for focus depths between 3 and 7 Mm is made, showing that 
approximately half of the observed signal originates within the first 2 Mm below 
the surface. A consequence of this is a rapid decrease (and reversal in some 
cases) of the travel-time perturbations with depth due to the contribution to 
the measurements of oppositely directed surface flows in neighboring convective 
cells. This confirms an earlier interpretation of similar effects reported from 
observations of supergranulation. 

Subject headings: Sun: helioseismology, interior 



1. Introduction 

The understanding of solar interior dynamics, as manifested in mass flows with a va- 
riety of temporal and spatial scales beneath the surface of the Sun, has been consider- 
ably advanced by both global- and local-helioseismic observations. Some recent reviews 
of th e progress in helioseismology o f flows (and other properties of the solar inte r ior) in- 



clude Christensen-Dalsgaardl (120021) . iThompson et all (120031 ) . iGizon & Birchl (120051 ). Miesch 



( 120051 ). and iThompsonl (120061 ) . In addition, progress has been made with increasingly 



improving numerical modeling procedures and computational resources. Examples of nu- 
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addition to providing physical insights into the origin and evolution of flows inferred through 
helioseismic observations and analysis, numerical models such as those cited above are being 
used to validate the observational and analysis tools themselves. The application of analy- 
sis methods to numerical (or "artificial") data can potentially l ead to improvemen t in both 
theory and observations of the phenomena being explored (e.g. IWerne et al.l 120041 ) . 



The mean solar rotation, and its variation with depth, latitude, and time, is detectable 



checks on the reliability of the analysis and modeling (e.g. 
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in local helioseismic methods is also obtained by comparing models, derived using differ- 



Giles et al. 1997. 1998: 1 
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Svanda et al. 


2007). Intercomparisons of methods and datasets can sometimes lead to the discovery of im- 



portant systematic effects or artifacts (IGonzalez Hernandez et al.l 120061 ) . 



The helioseismic study of "small-scale" structures within the Sun, defined here by wave- 
speed, flow, or magnetic perturbations in the solar interior with spatial scales less than a few 
tens of Mm is a particularly challenging pursuit. This is because the structures being probed 
have sizes which are of the order of, or perhaps only a few times larger than, the wavelength 
of the p- or /-modes used in the analysis. In this work we consider small-scale flows rep- 
resentative of solar supergranulation. The effects on mode frequencies and travel times of 
aco ustic waves caused by even smaller-scale flows (e.g. granu l ation) have also been explored 
(e.g iMurawski fc Roberts! Il993l ; iMurawski fc Goossenslll993l ; iPetrovay et al.l 120071 ). but are 
not considered here. Our simulations do not include magnetic fields and consequently we are 
unable to assess the effects of magnetic fields on the modeling or interpretation of helioseis- 
mic measurements of these flows. This important issue, however, continues to be extensively 
expl ored in the context of global helioseismology , asteroseismology, and local heliqseism ology 



(e.g. 
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Supergranulation has been a focus of local helioseismology for over a decad e 

Kosovichev fc Duvalj[l9~97l: buvall fc Gizon!l200ol:lGizon et al.lhoodbraun fc Lindsev 



20031 : Izhao & Kosovichevll2003l : iBraun et al.l 120041 ) . However, as noted by 



Gizon fc Birch! (120051 ) there is no definitive consensus on even the depth of the supergranula- 
tion phenomenon. Some inversions of p-mode travel times indicate the presence of a "return 
flow" (oppositely directed flows from what is o bserved at the solar surface) at various depths , 
but typically ~ 10 M m, below the photosphere (jDuvall et al.lll997l ; IZhao fc Kosovichevll2003l ). 
Braun et al.l (120041 ) have suggested that an observed change in sign in travel-time perturba- 
tions with increasing depth of the penetration of the modes may represent a contamination 
of the signal from neighboring supergranules. How this might effect the inversions, however, 
is not clear. I Zhao et al.l (120071 ) find that travel-time inversions on simulated supergranular- 
sized convection do not reproduce the simulation flow fields at depths below the photosphere 
greater than 3 Mm. How deep we are able to reliably measure subsurface flows due to super- 
granulation and other small-scale patterns remains an important, but unanswered, question. 

In this paper, we explore the prospects for helios eismic probing of small -scale flows in 
the solar interior by applying helioseismic holography (ILindsey fc Braunlll997l ) to recent nu- 
merical simulations of solar convection. To do this we compare the "observed" signatures 
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of flows (in this case, maps of the perturbations to p-mode travel times) obtained by ap- 
plying helioseismic holography to the surface of realistic solar simulations, to the expected 
("model") signatures obtained directly from the simulated flows. This comparison facili- 
tates an assessment of both the expected helioseismic signals and the uncertainties in the 
observations (i.e. departure of the observed signatures from the expected values). While the 
degree of similarity between observed and model signatures provides a critical validation of 
our specific helioseismic analysis our principle goal is more general than a test of a single 
method of observation or modeling procedure. In particular, the observational uncertainties 
estimated from this forward-modeling experiment are used to assess the general prospects 
for the seismic detection of small-scale flows in the solar interior. This is possible to the 
extent that 1) our simulation measurements have uncertainties characteristic of (or at least 
no worse than) local helioseismic observations in general, and 2) flows in the solar interior 
have properties similar to those in the simulations. The v alidity of the secon d condition is 
maximized by using realistic solar convection simulations ( Benson et al.l I2OO6I ) as described 
in § [2j To ensure the validity of the first of these conditions, we employ helioseismic holog- 
raphy in the "lateral vantage" , which (as explained below) is designed to utilize most of the 
p-modes propagating through a specific target point (focus) in the solar interior. We also 
compare the uncertainties estimated in the simulation observations with results obtained 
from the Michelson Doppler Imager (MDI) onboard the Solar and Heliospheric Observatory 
(SO HO). Our analysis and forward modeling procedures are described in more detail in §[3] 
and § H] respectively. The results of the comparison are shown in § [5], followed in § [6] by a 
discussion of the implications of these results. 



Simulations 



The 3D simulations of solar convection on s uperg ranular scales employed in this study 
were performed and described by iBenson et al.l (120061) and have already been used t o vali- 
date surface and subsurface flow diagnostics (jGeorgobiani et al.ll2007l ; IZhao et al.ll2007l ). The 
simulations are carried out in a box 48 by 48 Mm horizontally and 20 Mm deep. The conser- 
vation equations for density, momentum and internal energy are solved in conservative form, 
on a three-dimensional staggered mesh, using sixth order finite difference spatial derivatives 
and a low memory, third order Runge-Kutta time advance. The grid is uniform in horizontal 
directions and stretched in the vertical (stratified) direction. Horizontal boundary conditions 
are periodic, while top and bottom boundary conditions are open. Inflows at the bottom 
boundary have constant pressure, specified entropy and damped horizontal velocities. Out- 
flow boundary values are obtained by extrapolation. The code uses a tabular equation of 
state, that includes local thermodynamic equilibrium (LTE) ionization of the abundant el- 
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ements as well as hydrogen molecule formation, to obtain the pressure and temperature as 
a functio n of log density and int ernal energy per unit mass using the Uppsala atmosphere 
package (IGustafsson et al.lll975l ). Radiative heating/cooling is determined by solving the 
radiation transfer equation in both c ontinua and lines using the Feautrier met hod, assuming 
Local Thermodynamic Equilibrium (jNordlundlll982l ; IStein fc Nordlundll2003l ). 



The acoustic waves (p-modes) generated by convective motions in these simul ations 
have properties remarkably similar to observed solar p-modes (IGeorgobiani et al.ll2007l ). One 
notable exception is an enhancement of acoustic power at higher temporal frequencies relative 
to observed solar values. A fortunate consequence of this is the ability to perform and test 
helioseismic procedures with modes extending to higher frequencies than normally used with 
solar data. 

While power spectra of the convective motions in these simulations demonstrate the 
prese nce of structures over a broad range of spatial scales ( Benson et al.12006 ; Georgobiani et al 
20071 ). we expect that local helioseismic techniques are best suited for probing structures 
larger than the wavelength of the p-modes employed in the analysis. Figure [Th. shows the 
horizontal divergence of the 8.53 hr time-average of the simulated flows at a depth of 4 Mm 
below the surface. Figure [lb shows the same flow divergence after applying a Gaussian 
smoothing with a full- width-half-maximum (FWHM) of 4 Mm. This width is of the order of 
the wavelength of p-modes with a temporal frequency of 3 mHz just below the photosphere 
(or equivalently, the wavelength at 5mHz of modes at a depth of about 4 Mm). Evident in 
both panels (but especially in the smoothed panel) are cellular structures, the most promi- 
nent having diameters on the order of 20 Mm. Presumably these structures are the most 
amenable to helioseismic detection with mode wavelengths greater than 4 Mm, and serve 
as a proxy for solar supergranulation which has a somewhat larger scale (~ 30 Mm). The 
simulated flow structures are coherent to depths of about 12 Mm, below which appear weak 
counterflows (see Figure [2]). 

We use 8.53 hours of the simulations in our analysis, which is typical of solar helioseismic 
observations and is less than the lifetimes of individual solar supergranules (~ 25 hr). We 
apply helioseismic holography (see § E} to the vertical component of velocity sampled a 
height of 200 km above the continuum optical depth of unity in the background stratification 
(hereafter referred to as the "surface" of the simulations). We use flows sampled only every 
30 seconds in time, yielding 1024 total time samples. The original simulations are computed 
in a 500 by 500 horizontal grid (giving a horizontal resolution of 96 km) and are for our 
purposes degraded in resolution by a factor of two in each direction by pixel averaging. 
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3. Analysis 



Helioseismic holography (hereafter HH) is a method based on the phase-coherent imag- 
ing of the solar interior aco ustic field. It computationally extrapolates the surface acoustic 
field into the solar interior (ILindsey fc Braunl 119971 . 120001 ) to estimate the amplitudes of the 
waves propagating into and out of a focus point at a chosen depth and position in the solar in- 
terior. These amplitudes, called the ingression and egression, are es timated by a convolutio n 
of the surface oscillation signal with appropriate Green's functions (ILindsey fc Braunl 120001 ) . 
Here, HH is performed in the wavenumber-frequency (Fourier) domain using a set of fixed 
frequency bandpasses with a width of 1 mHz and centered at frequencies, u , of 3, 4, 5, and 
6 mHz. While previous applications of HH to study solar flows have been limited, by signal 



to-no i se considerations, to frequencies equal to and below 5.5 mHz (e.g. iBraun fc Lindsey 
20031 ; IBraun et al.l 12004 ). we are motivated by the presence of high frequency waves in the 
simulations to include a 6mHz bandpass in our analysis here. 



The HH analysis is performed in the wave-mechanical formulation (ILindsey fc Braun 



20041 ) . The method employed for horizontal flow diagnostics is base d on the egressions and 
ingressions computed in the lateral vantage (ILindsey fc Braunl 120041 ) employing pupils span- 
ning 4 quadrants extend ing in different directions (east, west, north and south) from the 
focus (IBraun et al.ll2004l ). In the lateral vantage, the p-modes sampled by the pupil prop- 
agate through the focal point in directions inclined up to ±45° from the direction parallel 
to the surface (Figure [3]). The antisymmetric phase shift, the difference in the phase per- 
turbation of waves traveling from one pupil to its opposite and the phase perturbation of 
waves traveling in the reverse direction, is sensitive to horizontal flows near the focus. In 
general, the phase perturbation S(p A is related to an equivalent travel-time perturbation by 
St = 8<j) /2TT1/Q. The travel time perturbations, St x and 5r y derived from antisymmetric 
phase shifts computed between the east/west and north/south quadrant pairs, respectively, 
provide the HH signatures sensitive to the two components of the horizontal flow. The sign 
of the travel-time perturbations are such that a positive velocity component will produce a 
negative value of St (i.e. a reduction in the mode propagation time). The lateral- vantage 
geometry samples more than 70% of the wave modes which pass through the focus. The 
remaining waves, propagating more vertically than the waves appearing in the pupil, are 
substantially less sensitive to horizontal flows. 

Table [T] lists the focus depths and the pupil radii used in lateral- vantage HH. The pupil 
radii are defined from ray theory. The range of (spherical-harmonic) mode degrees (£) at 
4 mHz, selected by each pupil, is listed in the table. The lower £ value denotes the modes 
propagating at ±45° from the horizontal direction which propagate through the focus and 
reach the surface at either the inner or outer pupil radius. The highest £ value listed indicates 
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modes propagating horizontally through the focus. The mode degrees selected by the other 
frequency bandpasses scale approximately with u . Dispersion effects predict slight variations 
in the pupil radii with wave frequency. However, experimentation shows that the measured 
travel-times do not vary with similar adjustments of the pupil radii sufficiently to affect the 
comparisons and conclusions discussed here. Consequently, for convenience and simplicity, 
we use the same pupil for all frequency bandpasses for a given depth focus. 

Comparisons of travel-time maps made with and without the use of phase-speed filters 
show that these filters are of marginal utility when applied to data subsequently sampled over 
wide pupils of the sort employed in lateral- v antage HH. We not e that this is not the case for 



much of time-distance helioseismologv (e.g. iDuvall et al.l 119971 ; IZhao. Kosovichev fc Duvall 
200ll ; IZhao et al.l 120071 ) or for HH performed with narrow pupils ( IBraun fc Birch 2006), 



where the use of phase-speed filters appears to provide a dramatic improvement in the 
measurements. While a reduction in scatter is sometimes evident with the use of these filters 
in lateral vantage HH, we find that, for both solar data and the present simulated data, the 
effect is comparable to what can be achieved by spatially smearing the unfiltered travel-time 
perturbation maps. In fact, an understanding of the consequences of smearing the travel- 
time maps is critical for interpreting our comparisons between observed and modeled results. 
With these considerations in mind, our strategy is to use observed travel-time maps made 
without phase-speed filters, and explore explicitly the direct consequences to our comparisons 
of spatial averaging of the unfiltered maps by varying amounts. 

The focus depths chosen to perform the "observations" extend down to about 8 Mm 
below the simulation surface. At focus depths larger than this, egression-ingression corre- 
lations cease to exhibit meaningful results, even though reliable correlations are routinely 
observed at these depths in the Sun. This failure may be due to the lower boundary of the 
simulations at 20 Mm. Even though the lower boundary of the simulation is twice the depth 
of the focus depth where problems appear to arise, we note that the lateral-vantage geometry 
samples p-modes which penetrate a considerable depth below the focus, and any reflection 
or absorption of these waves at the lower boundary will likely adversely affect our observa- 
tions. The horizontal periodic boundaries used in the simulations ensure that HH travel-time 
measurements, performed in the Fourier domain, are usable over the full horizontal extent 
of the simulations. Thus, the data is "infinite" (but periodic) horizontally, and travel-time 
perturbation maps are free from any detrimental effects caused by an abrupt termination of 
data at the edges of the domain, such as typically experienced with helioseismic observations 
of the Sun. 
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Forward Models 



In this section we use the Born approximation approach of iGizon fc Birch! (120021 ) to 
estimate the HH signatures that would be expected from the flows in the simulations. The 
range of validity of the Born approximation for three-dimensional time-dependent flows in 
stratified models is not known. 



Birch fc Felderl (120041 ) studied the validity of the Born approximation as applied to time- 
distance helioseismology of simple steady flows in 2d unstratified models. In these models, the 
Born approximation is typically seen to be valid when travel-time shifts are a small fraction 
of the wave period. The validity of the Born approximation for scattering from steady sound- 
speed inhomogeneiti 
helioseismology (e.g 



Fan et al. 


1995: 


Hune et al. 


2000: 


Birch et al. 


2001; 


Baie et al. 


2003) 



In these studies, it is again typical that the Born approximation is valid when travel-time 
shifts are small compared to the mean wave period (though it can sometimes be valid well 
past this limit). The extension of this general conclusion to the current work, which involves 
strong time-dependent small-scale flows, is not at all certain. 

Comparison between the Born approximation estimates described in this section and 
the actual travel times (§ \5§ will provide an important estimate of the validity of the Born 
approximation for complicated solar-like flows. In particular, we will be able to study the 
extent to which the HH signatures can be predicted using only the time average of the flow. 



As described in detail by lGizon &: Birch! (120021 ) the computation of the linear sensitivities 
of local helioseismic measurements requires, in general, two steps. The first step is the 
computation of the linear sensitivity of the measurement procedure (e.g. the measurement 
of travel times) to small changes in covariance of the wavefield. The second step is the 
computation of the change in the wavefield covariance introduced by small changes in the 
model of the solar interior (e.g. the int roduction of flows or sound-speed variations). This 
general procedure has been a pplied by IBirch fc Gizonl (120071 ) to the case of time-distance 
measurements of flows and by IBirch et al.l (120071 ) to the case of ring-diagram measurements 
of flows. 



The application of the IGizon fc Birch! (120021 ) recipe to the case of HH is relatively 
straightforward. The result is a set of linear sensitivity functions (kernels), K, which give 
the linear sensitivity of the anti-symmetric phase 4> A to small amplitude steady flows v(x), 



iA 



dx K(x) • vfx) 



(1) 



where x denotes three-dimensional position in the model and the integral is taken over the 
entire three-dimensional domain of the simulation. In the computation of the kernel functions 



- 9- 



we use the norma l-mode Greens functions, source model, and dampi ng model described by 
Birch et al.l (12004 ) . The background stratification is given by Model S (jChristensen-Dalsgaard et al 
1996 



One limitation of the normal-mode Green's functions of iBirch et al.l (120041 ) is that they 
employ the upper boundary cond ition that Lagrangi an pressure perturbation vanishes at the 
top of Model S. As discussed by IBirch et al.l (120041 ). this boundary condition is reasonable 
for waves at frequencies well below the acoustic cutoff frequency. For frequencies near the 
acoustic cutoff frequency, this boundary condition is likely not appropriate. As a result, the 
kernels described in this section are not intended to model HH measurements at or above 
five mHz. 



5. Results 

5.1. Frequency and Depth Variations of Observed Travel-Time Perturbations 

An examination of the observed travel-time perturbation maps reveals a fair degree of 
correlation amongst different frequency bandpasses for a given focus depth, particularly for 
the shallower depths (< 5 Mm) selected for analysis. An example, at the 0.7 Mm focus 
depth, is shown in Figure HI Correlation coefficients between perturbation maps at different 
frequency were computed after applying Gaussian smoothing with varying FWHM to the 
maps. As expected, the correlations improve with increased smoothing. With a Gaussian 
FWHM of 4 Mm, the correlation coefficients are of the order of 0.7 for the shallowest depths 
and decrease to about 0.5 at a depth of 6 Mm. At depths of 7 and 8 Mm, the 1 mHz 
bandpass maps are essentially uncorrelated with each other. As with analogous travel-time 
observations in the Sun, the measurement noise from these simulations is most likely due 
to realization noise caused by stochastic excitation of the p-modes. In the solar case, the 



prop erties of the measurement noise have been extensively explored (e.g. iGizon fc Birch 



20041 ). 



At depths where comparisons between frequencies are meaningful, we do not detect any 
clear systematic variation of the travel-time perturbations with frequency. Typically, the 
slope of a linear fit of the travel times between any two frequencies differs from unity by 
less than 10%. A somewhat surprising result is that the variance of the 1 mHz bandpass 
maps, defined as the average over all pixels of the squared difference between each map and 
an average of maps over all four bandpasses (hereafter referred to as a frequency-averaged 
map), is essentially the same for all frequencies between 3 and 6 mHz. 



As expected, there is a good correspondence, particularly at shallower depths, between 
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the travel-time maps and the relevant component of the horizontal velocity. This is partic- 
ularly evident in comparisons between the frequency-averaged perturbation maps and the 
simulation flows (e.g. compare Figure @£ and Figure Hf). However, it is also apparent that 
the ratio of the magnitude of the travel-time perturbations to their values at the shallowest 
focus depth (shown by the dotted lines in Figure [2]) falls off considerably faster with focus 
depth than a similar ratio of the simulation flow magnitudes (solid lines). 



5.2. Comparisons with Model Travel-time Perturbations 

Using the methods described in §@]we derive the model travel-time perturbations (5r mod ) 
for comparisons with the observed values. Because of the restriction to low frequency modes 
described in § |4]we compute <5r mo d only for the lowest frequency bandpass {uq = 3 mHz). 
Figure [5] shows a comparison between the frequency-averaged travel-time perturbations and 
the 3-mHz model for three focus depths; 3, 5, and 7 Mm below the surface. Also shown are 
the residuals after subtracting the model perturbations from the observations. There is good 
agreement between the model and observations for the shallower depths. At greater depths 
the similarity becomes less apparent. The patchy appearance of the residual maps is likely 
due to the properties of t he realization noise w hich is coherent over distances of about half 
of a p-mode wavelength (IGizon fc Birch! 12004 ) . The similarity of the residuals at different 



focus depths follows from the highly overlapping set of modes used in the analysis for the 
different depths (see Tabled]). We define the square of the deviation from the model, cx^od; 
of each observed travel-time perturbation map (5r b s ) as: 



a mod = 5Z^ Tobs ~ 5r mod) 2 , (2) 

where the summation is over all N pixels in the maps. 

To see how a mod depends on smoothing, we spatially convolve each observed perturba- 
tion map by two-dimensional Gaussian functions with varying FWHM. Figure [6] shows the 
resulting <7 mo d as a function of FWHM. The different colors indicate different focus depths, 
and the results are shown for both the 3 mHz bandpass (dotted lines) and the average over 
the four frequency bandpasses (solid lines). The deviation in the frequency-averaged pertur- 
bation maps is reduced from the single 3 mHz bandpass by a factor of two (i.e. the decrease 
from the dotted to solid lines) as expected from independent measurements with similar 
uncertainties (although this improvement declines with more smoothing). 

If each pixel represented an independent sample of the travel-time perturbation with 
an uncertainty which is normally distributed, we would expect cr mo d to decrease with the 
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inverse of the FWHM. Figure [6] shows, however, a considerably weaker decrease of the 
deviation, most likely due to the finite wavelengths of the modes used to produce the maps 



( IGizon Sz Birchll2004j ) . In other words, the maps are already "smoothed" by finite- wavelength 
effects which can be readily seen in the raw perturbation maps. With values of FWHM 
more than about 10 Mm, the effects of smoothing increase the departure from the model 
signatures, since the level of smoothing compromises the ability to resolve the dominant flow 
structures (e.g. Figured]). When "over-smoothed" in this way, there is less (and sometimes 
no) improvement achieved by frequency averaging. The values of cr mod are generally similar 
for all of the depths shown, although the deeper two measurements can apparently withstand 
slightly greater smoothing to obtain smaller deviations from the model. 

We have examined the statistics of the travel-time measurements in the independent one 
mHz-wide frequency bandpasses by computing the standard deviation of the measurements 
in each pixel from the mean over the four frequency bandpasses. Some caution is warranted 
in directly comparing the pixel- average of this quantity, which we denote a u , with <T mo d when 
the data is smoothed. For example, substantial smoothing can reduce the difference between 
frequency bandpasses to an arbitrarily small value, even while both differ substantially from 
the model. A correspondence between these two quantities is expected however, in the 
limit of little or no smoothing, if there are no systematic differences between frequency 
bandpasses of the perturbations or their uncertainties. The values of a v from the maps 
of <5r bs, smoothed with a Gaussian FWHM only minimally larger than the pixel size, are 
indicated by the diamonds in Fig. The correspondence is particularly striking for depths 
of 3 and 5 Mm. At a depth of 7 Mm, the ratio a u /a mo ^ is about 0.8. 



5.3. Comparisons with MDI Observations 

It is of interest to compare the properties of the simulation observations with results 
obtained with actual solar observations. An 8.53 hour time series, with 60 second cadence, 
of solar images was processed identically to the simulated data reported here. The data 
set consists of 512 Dopplergrams obtained by SOHO/MDl JScherrer et. al.l Il995h on 1999 



June 1. The MDI Dopplergrams were interpolated onto Postel's projections, tracked with 
a Carrington rotation rate and centered on longitude L = 306° and latitude B = 7°. The 
projected frames consist of a grid of 512 by 512 pixels with a resolution of 828 km per pixel, 
obtained after a two by two pixel averaging of the original high resolution pixels. Travel- 
time perturbations were computed, with the same HH analysis as the simulations, for focus 
depths of 3, 5, and 7 Mm below the surface over 1-mHz temporal bandpasses centered at 
3, 4, and 5 mHz. In general, HH travel-time measurements from solar data, unlike the 
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simulations, are not reliable at higher frequencies. Figure [7] shows cropped (48 by 48 Mm) 
travel-time perturbation maps, averaged over all frequency bandpasses, of a region close to 
the center of the solar disk, compared with maps obtained from the simulated data at the 
same focus depths. It would be difficult to identify, from general appearances alone, which 
measurements came from the artificial data and which came from solar data (except perhaps 
for the finer pixel size noticeable in the simulations). 

To assess the noise in the MDI observations, we compute the pixel- averaged standard 
deviation across the three frequency bandpasses. For these measurements, only a minimal 
smoothing (FWHM w 1.2 Mm) is applied. The results are shown by the X symbols in 
Figure [6j These values agree very well with the single-bandpass model-deviations of the 
simulation measurements after smoothing with the same Gaussian function. We conclude 
that the magnitude of the uncertainties in HH travel-time perturbation measurements made 
with 8 hr of solar observations are very similar to those obtained with the simulated data 
analyzed here. 



5.4. Depth Contribution of Signal 



We use the forward models to examine the relative contribution with depth to the 
expected travel-time perturbations from the flow structures. To successfully infer the mag- 
nitude of subsurface flows from the travel-time observations, we will very likely require that 
the expected perturbations due to flows at a desired focus depth exceed the noise present in 
the observations. Thus we need to extract the relevant portion of the expected signatures (i.e. 
due to the flow near the target focus) from the total travel-time perturbations (which result 
from flows over the entire depth range sampled by the modes used in the measurements). 

Assessing the contribution to the travel-time signatures of flows near the surface are 
of particular interest. This is motivated by the desire to understanding the rapid decrease 
in the observed travel-time perturbations with focus depth as exhibited in Figure [2J Pre- 
vious HH analyses and modeling of supergranulation observed with SOHO/MDl have indi- 
cated the possibil ity of a rapid decrease, and eventual rever sal, of the travel-time perturba- 
tions with depth (iBraun fc Lindseyll2003l ; iBraun et al.ll2004l ). An explanation was offered by 
Braun fc Lindseyl (120031 ) that the results were consistent with an increase with focus depth 
of the near-surface contribution to the travel-time perturbations from oppositely directed 
flows in neighboring supergranule cells. These near-surface contributions increased with fo- 
cus depth because the pupils needed for lateral-vantage HH increase in radii with depth (e.g. 
Table [I]) eventually exceeding the size of the supergranular cells. Forward models a ssuming 
simple, shallow, cellular flows were able to reproduce this effect (IBraun et al.l 120041 ). and it 



is of considera ble interest to see if such behavior exists for the flow patterns present in the 
simulations of iBenson et al.l (120061 ). 



To assess the depth-dependence of the contribution to the travel-times due to subsurface 
flows we perform a horizontal convolution of the kernels (§ HJ) with the x component of the 
flow fields. Figure [8] shows a vertical slice of the x component of the flows, a slice through the 
kernel functions (for a focus depth of 3 Mm), and a slice through the horizontal convolution 
of the two. Figs M and [TU] show the the same for focus depths of 5 and 7 Mm, respectively. 
What is readily apparent in these figures is that, while the flows fall off gradually with depth 
(Figure EJa,), the contribution to the travel-time perturbations (Figs [Hb, [9_b, and [TUb) due to 
near-surface layers (with depths from 0-1 Mm) is disproportionately large. This results 
from the strong near-surface sensitivity of the travel-time kernels (Figs [8b, [Hb, and [TUb). It 
is noteworthy that the root-mean-square (RMS) of the travel-time perturbation integrated 
from the surface to a depth of 2 Mm is, for all three depths shown, roughly the same as the 
RMS of the perturbation integrated from 2 Mm to the bottom of the simulation, where the 
RMS is computed over the full horizontal domain. 

For what follows we define the "target signal" to be the contribution to the travel- 
time perturbation over a "target-range" of depths about a given focus depth. This range is 
intended to be proxy for the depth resolution of any modeling of the observations (through 
either forward or inverse methods). For a depth of 3 Mm below the surface, the integration 
is carried out from 2 Mm to 6 Mm. For the focus of 5 Mm, the target range is from 4 Mm 
to 8 Mm and for 7 Mm, the target range used is 5 to 10 Mm. These ranges were selected 
by a visual inspection of the sensitivity functions (e.g. FigsEb", [9Jd, and [TUb). Figs [Hi, [9ji, 
and [TUB, show plots of the total travel-time perturbations (blue lines), the target signal (red 
lines), and the contribution away from the target (black lines), over a slice of the data cube. 
Note that Figs. [UJ and (TUJ indicate a positive contribution to the travel-time perturbations 
at x ~ 30 Mm from the near surface, which is opposite in sign to what is expected for the 
rightward-directed flow at this location (Figure [8b). This contribution is due to flows located 
in the near surface lobes of the kernel which have opposite sign from the flow at the focus. 
These measurements offer support to the suggestion that the rapid decrease (and potential 
reversal in sign) of the total travel-time perturbation with increasing depth results from the 
near-surface contribution of nearby (and oppositely directed) flows. 

To estimate the signal-to-noise ratio (SNR) required for the detection and modeling of 
the subsurface flows in the simulations we use the RMS of the target signals as defined above. 
For the noise levels, we use the optimal (i.e. minimum) values of cr mo d from Figure [UJ for 
each depth. For 3 Mm, this is 5 seconds, and for 5 and 7 Mm, this is 4.5 seconds. Table [2] 
shows the values of the RMS for the total travel-time perturbation (column 2), the RMS of 
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the target signal (column 3) and the resulting SNR (column 4) for the three depths listed 
in column 1. These results do not substantially change if the "target" is adjusted to include 
flows at deeper layers. For example, extending the lower depth limit of each target range to 
the bottom of the simulation box increases the RMS target signal and the resulting SNR by 
approximately 5-15% for the depths listed in Table [2j It is clear that the SNR is too small 
for detecting subsurface flows in these simulations at and below 5 Mm with 8 hr of data. 
Note that a SNR of one implies the observed travel times are consistent, within the noise, 
with no flows. A SNR greater than unity is required for an unambiguous detection of a flow. 



6. Discussion 

A major finding of this work is that lateral-vantage holography recovers the model 
travel times about as well as one should expect. In other words, the RMS of the residuals 
(observed minus model) is very close to the noise levels as deduced by the variance between 
independent frequency bandpasses. We can also infer from the similarity of results across 
the different bandpasses that the validation of the method applies across a wide range of 
temporal frequencies. In addition, we have gained a fair degree of confidence in the use of 
the Born approximation for modeling p-mode travel-times for flows similar to the the ones 
simulated here. We note, however, that models have only been computed in one frequency 
bandpass and that further effort is needed to extend the modeling to higher frequencies. One 
important issue not addressed here is the potential improvement in the spatial resolution 
of flows (or other perturbations) deduced with observations made with increasing temporal 
frequency of the p-modes employed in the analysis. This is plausible if the resolution scales 
with the mode wavelength which, for a mode propagating through a fixed focus, varies 
inversely with the temporal frequency. 

We have made a quantitative assessment of the contribution of near-surface layers to 
travel-time measurements as a whole. A major consequence of our result is a plausible ex- 
planation for the relatively rapid decrease with depth (and eventual sign reversal) of the 
travel time perturbations, which has also been detected in MDI observations of solar super- 
granulation. This trend is not necessarily detrimental to successful modeling of subsurface 
flows, given a sufficient SNR in the travel time measurement and the assumption that the 
near-surface perturbations are well understood. Most modeling efforts typically include an 
assessment of flows at a range of depths, including the near-surface regions. In the simu- 
lations performed here, the horizontal flows decrease gradually with depth (Figs [2] and [5k>). 
However, strong shear layers immediately below the photosphere or the existence of surface 
perturbations due to magnetic fields (neither of which is present in these simulations) could 
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enhance the near-surface contributions and potentially complicate the modeling of flows. 
Some possible complications include strong surface perturbations which are not easily mod- 
eled under assumption of the Born approximation, or that might req uire the use of "surfac e 



terms" which are not routinely employed in local helioseismology (see lBraun &: Birch! 120061 ). 



Numerical simulations provide an excellent means of testing these possibilities. 

In as far as these simulations predict the types of flows likely to be present in Sun, it is 
appropriate to extrapolate our results to evaluate the prospects for the detection and mod- 
eling of actual solar flows. We assume similar noise levels, as confirmed by the comparisons 
performed in § 15.31 An increase in the temporal duration of the observations from 8 to 24 
hrs of observations, or roughly the lifetime of the super granules, should produce an increase 
in the SNR by a factor of y/3. On the other hand, our experience has shown that only 
the 3, 4, and 5 mHz frequency bandpasses can be used successfully in solar measurements, 
so that the nominal increase in SNR is a more modest factor of about 1.5. The resulting 
SNR values are shown in column 5 of Table El The general result, that supergranule-sized 
flows are essentially undetectable using current methods below depths around 5 Mm, is not 
substantially changed. 

It is fair to point out that our pessimistic conclusions about detecting small-scale sub- 
surface flows below depths of 5 Mm are based on idealized conditions and assumptions. To 
the extent that these conditions may not be representative of actual or proposed solar mea- 
surements and modeling it is plausible that our conclusions may not be fully relevant. On 
the other hand it has been our goal to err, if at all, on the side of optimism in computing 
the SNR values quoted here. Our noise estimates are inferred from forward models con- 
structed directly from the subsurface flows present in the simulations, and thus represent an 
idealized, "best-case," scenario. Thus, the SNR values quoted do not address possible errors 
introduced by mo deling limitations or systematic effects encountered in solar measurements. 



Zhao et al.l (120071 ) have performed inversions of time-distance measurements performed on 
this simulation and find only weak correlations between the actual and inverted flows at 
depths 3 Mm and greater below the surface. This appears to be consistent with our findings 
regarding SNR estimates, although the methods employed are different. 

While the use of lateral- vantage HH is a fairly efficient means of measuring the effects of 
the horizontal component of the flow field (see §[3]), we recognize that additional information 
may be used to infer the properties of subsurface flows. These include using helioseismic 
measurements of vertical flows and the assumption of mass conservation as an additional 
constraint on horizontal flow components, or matching subsurface flows to those observed 
at the surface through other means (e.g. direct Doppler measurements or feature-tracking). 
The degree to which this additional information can be used to substantially improve our 
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ability to measure flows deeper than a few Mm below the photosphere, however, remains to 
be determined. We also recognize that solar supergranulation has a spatial scale somewhat 
larger than the cells explored in these simulations, which may permit somewhat more spatial 
smearing and subsequent reduction in noise. Some improvement might be expected by 
increasing the temporal duration of the observations, since some supergranules last as long 
as several days. And obviously, many subsurface dynamic phenomena, both known and 
unknown, will still be accessible after temporal and or spatial averaging. It is also important 
to keep in mind the spatial and temporal scales of the flows for which our findings are relevant. 
Clearly, SNR values far greater than unity are possible, even at substantial depths below the 
photosphere, for measurements of large-scale or long-lived structures such as meridional and 
zonal flows. 



Birch et al.l (120061 ) presented a method to combine travel-time measurements, from 
time-distance helioseismology, of thousands of supergranules, allowing the exploration of 
supergranule-averaged flows. In general, a parameterization of small-scale flows (in terms 
of statistical moments or Fourier components, for example) m ight be more amenab le to 



modeling than the reconstruction of individual flow vectors (e.g. I Woodardll2006l . 120071 ) 



Forward modeling efforts with both solar and artificial data appear to offer considerable 
utility for understanding both advantages and limitations of helioseismic methods. It would 
be useful to extend the analysis presented here to address a variety of additional issues. 
Understanding the trade-off between maximizing spatial resolution in modeling perturbations 
and minimizing the uncertainties in the models is a major issue which can be explored 
through forward modeling. For example, it should be straightforward to measure the signal 
and noise values for different spectral components (in the spatial domain) of flows and other 
perturbations to address this. 

Finally, we emphasize the importance of numerical simulations for helioseismic analysis. 
Some improvements are now becoming available or are under preparation. These include 
simulations which cover larger areas and extend deeper into the solar interior and are carried 
out in spherical geometry. There is also a need to compliment realistic simulations (which 
give the best current estimate of what's likely to be found in the solar interior) with the 
ability to prescribe and adjust perturbations (in order to fully test observing and modeling 
methods and to allow for surprises in the Sun). The addition of magnetic fields to numerical 
simulations will be of particular utility to test both their influence on convective and other 
flows, and their influence on local helioseismic measurements. 



DCB and ACB are supported by funding through NASA contracts NNH05CC76C and 
NNH04CC05C, NSF grant AST-0406225, and a subcontract through the HMI project at 



-17- 



Stanford University awarded to NWRA. The numerical simulations were carried out with 
support by NASA grant NNG04GB92G and NSF grant AST-0605738. 

REFERENCES 

Baig, A. M., Dahlen, F. A., k Hung, S.-H. 2003, Geophysical Journal International, 153, 
467 

Basu, S., k Antia, H. M. 2000, Sol. Phys., 192, 469 

Basu, S., Antia, H. M., & Tripathy, S. C. 1999, ApJ, 512, 458 

Benson, D., Stein, R., k Nordlund, A. 2006, Astronomical Society of the Pacific Conference 
Series, 354, 92 

Birch, A., Duvall, T. L., Gizon, L., k Jackiewicz, J. 2006, AAS/Solar Physics Division 
Meeting, 37, #05.05 

Birch, A. C, k Felder, G. 2004, ApJ, 616, 1261 

Birch, A. C, & Gizon, L. 2007, Astronomische Nachrichten, 328, 228 

Birch, A. C, Gizon, L., Hindman, B. W., and Haber, D. A. 2007, ApJ, submitted 

Birch, A. C, Kosovichev, A. G., k Duvall, T. L., Jr. 2004, ApJ, 608, 580 

Birch, A. C, Kosovichev, A. G., Price, G. H., k Schlottmann, R. B. 2001, ApJ, 561, L229 

Bogdan, T. J. 2000, Sol. Phys., 192, 373 

Braun, D. C, k Birch, A. C. 2006, ApJ, 647, L187 

Braun, D. C, Birch, A. C, k Lindsey, C. 2004, ESA SP-559: SOHO 14 Helio- and Astero- 
seismology: Towards a Golden Future, 14, 337 

Braun, D. C, k Fan, Y. 1998, ApJ, 508, L105 

Braun, D. C, k Lindsey, C. 2003, ESA SP-517: GONG+ 2002. Local and Global Helioseis- 
mology: the Present and Future, 12, 15 

Chou, D.-Y., k Ladenkov, O. 2005, ApJ, 630, 1206 

Christensen-Dalsgaard, J., et al. 1996, Science, 272, 1286 



Christensen-Dalsgaard, J. 2002, Reviews of Modern Physics, 74, 1073 
Duvall, T. L., Jr., et al. 1997, Sol. Phys., 170, 63 
Duvall, T. L., Jr., & Gizon, L. 2000, Sol. Phys., 192, 177 
Fan, Y., Braun, D. C, & Chou, D.-Y. 1995, ApJ, 451, 877 

Georgobiani, D., Zhao, J., Kosovichev, A. G., Benson, D., Stein, R. F., & Nordlund, A. 
2007, ApJ, 657, 1157 

Giles, P. M., Duvall, T. L., Jr., & Scherrer, P. H. 1998, Structure and Dynamics of the 
Interior of the Sun and Sun-like Stars, 418, 775 

Giles, P. M., Duvall, T. L., Jr., Scherrer, P. H., & Bogart, R. S. 1997, Nature, 390, 52 

Gizon, L., & Birch, A. C. 2002, ApJ, 571, 966 

Gizon, L., & Birch, A. C. 2004, ApJ, 614, 472 

Gizon, L., & Birch, A. C. 2005, Living Reviews in Solar Physics, 2, 6 

Gizon, L., Duvall, T. L., Jr., & Larsen, R. M. 2000, Journal of Astrophysics and Astronomy, 
21, 339 

Gizon, L., Duvall, T. L., & Schou, J. 2003, Nature, 421, 43 

Gonzalez Hernandez, I., Komm, R., Hill, F., Howe, R., Corbard, T., & Haber, D. A. 2006, 
ApJ, 638, 576 

Gustafsson, B., Bell, R. A., Eriksson, K., & Nordlund, A. 1975, A&A, 42, 407 

Haber, D. A., Hindman, B. W., Toomre, J., Bogart, R. S., Thompson, M. J., & Hill, F. 
2000, Sol. Phys., 192, 335 

Haber, D. A., Hindman, B. W., Toomre, J., Bogart, R. S., Larsen, R. M., & Hill, F. 2002, 
ApJ, 570, 855 

Hanasoge, S. M., et al. 2006, ApJ, 648, 1268 

Hughes, S. J., & Thompson, M. J. 2003, GONG+ 2002. Local and Global Helioseismology: 
the Present and Future, 517, 307 

Hung, S.-H., Dahlen, F. A., & Nolet, G. 2000, Geophysical Journal International, 141, 175 



-19- 



Jensen, J. M., Olsen, K. B., Duvall, T. L., Jr., & Jacobsen, B. H. 2003, ESA SP-517: 
GONG+ 2002. Local and Global Helioseismology: the Present and Future, 12, 319 

Khomenko, E., & Collados, M. 2006, ApJ, 653, 739 

Kosovichev, A. G., & Duvall, T. L., Jr. 1997, SCORe'96 : Solar Convection and Oscillations 
and their Relationship, 225, 241 

Lindsey, C. & Braun, D. C. 1997, ApJ, 485, 895 

Lindsey, C. & Braun, D. C. 2000, Sol. Phys., 192, 261 

Lindsey, C. & Braun, D. C. 2004, ApJS, 155, 209 

Mansour, N. N., Kosovichev, A. G., Georgobiani, D., Wray, A., & Miesch, M. 2004, ESA 
Special Publication, 559, 164 

Miesch, M. S. 2005, Living Reviews in Solar Physics, 2, 1 

Murawski, K., & Goossens, M. 1993, A&A, 279, 225 

Murawski, K., & Roberts, B. 1993, A&A, 272, 601 

Nordlund, A. 1982, A&A, 107, 1 

Parchevsky, K. V., & Kosovichev, A. G. 2007, ArXiv Astrophysics e-prints, 
[arXrv:astro-ph/0612364| 

Petrovay, K., Erdelyi, R., & Thompson, M. J. 2007, Sol. Phys., 240, 197 

Scherrer, P. H. et al. 1995, Sol. Phys., 162, 129 

Shelyag, S., Erdelyi, R., & Thompson, M. J. 2006, ApJ, 651, 576 

Shelyag, S., Erdelyi, R., & Thompson, M. J. 2007, A&A, 469, 1101 

Stein, R. F., & Nordlund, A. 2000, Sol. Phys., 192, 91 

Stein, R. F., & Nordlund, A. 2003, Stellar Atmosphere Modeling, 288, 519 

Svanda, M., Zhao, J., & Kosovichev, A. G. 2007, Sol. Phys., 241, 27 

Thompson, M. J. 2006, Royal Society of London Philosophical Transactions Series A, 364, 
297 



-20- 



Thompson, M. J., Christensen-Dalsgaard, J., Miesch, M. S., & Toomre, J. 2003, ARA&A, 
41, 599 

Tong, C. H., Thompson, M. J., Warner, M. R., & Pain, C. C. 2003, ApJ, 593, 1242 

Werne, J., Birch, A., & Julien, K. 2004, ESA SP-559: SOHO 14 Helio- and Asteroseismology: 
Towards a Golden Future, 14, 172 

Woodard, M.F., 2006, ESA SP-624: Proceedings of SOHO 18/GONG 2006/HELAS I, Be- 
yond the spherical Sun, 18 

Woodard, M.F., 2007, ApJ, submitted 

Zhao, J., Georgobiani, D., Kosovichev, A. G., Benson, D., Stein, R. F., & Nordlund, A. 
2007, ApJ, 659, 848 

Zhao, J., & Kosovichev, A. G. 2003, GONG+ 2002. Local and Global Helioseismology: the 
Present and Future, 517, 417 

Zhao, J., k Kosovichev, A. G. 2004, ApJ, 603, 776 

Zhao, J. & Kosovichev, A. G. & Duvall, T. L., Jr. 2001 ApJ, 557, 384 



This preprint was prepared with the AAS IATgX macros v5.2. 



-21 - 



Table 1. Pupil radii. 



Depth 


Pupil radii 


i @ 


4mHz 


(Mm) 


(Mm) 






0.7 


0.7 - 11.1 


980 • 


- 1370 


1.5 


1.2 - 12.5 


850 • 


- 1190 


2.3 


1.8 - 13.9 


750 ■ 


- 1060 


3.0 


2.1 - 14.6 


680 


- 960 


4.0 


2.8 - 16.7 


590 


- 840 


5.0 


3.5 - 18.8 


530 


- 740 


6.0 


3.5 - 24.4 


470 


- 670 


7.0 


4.2 - 29.9 


420 


- 600 


8.3 


4.9 - 39.0 


370 


- 520 
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Table 2. Signal-to-Noise Estimates 



Depth 


RMS total signal 


RMS target 


SNR 


SNR 


(Mm) 


(s) 


(s) 


8 hr sims. 


24 hr Sun 


3.0 


15.6 


8.7 


1.7 


2.6 


5.0 


10.5 


3.6 


0.8 


1.2 


7.0 


7.4 


1.7 


0.4 


0.6 
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Fig. I. — a) The horizontal divergence of the time average of the simulated flows at a 
depth of 4 Mm below the surface. Positive values of the divergence, identified by bright 
regions, correspond to outflows, b) the same flow divergence smeared with a two dimensional 
Gaussian function with a FWHM of 4 Mm. 
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Fig. 2. — The ratio with respect to near-surface values of the simulated horizontal flows (solid 
lines) and the measured holographic travel-time perturbations (dotted lines) as a function of 
depth. The blue and red lines indicate the measurements for the x and y components of each 
quantity respectively. For the travel time perturbations, the depth corresponds to the depth 
of the focus in the lateral vantage. For both flows and travel time perturbations, the ratio 
is defined to be the slope of a linear least-squares fit between the relevant quantity at the 
indicated depth with the quantity evaluated at a "near-surface" depth. For the simulated 
flows this depth is the surface of the simulations. For the travel-time perturbations this 
depth is 0.7 Mm below the surface. 
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Fig. 3. — Lateral-vantage holography with the focus placed at a depth of 7 Mm below 
the surface. Solid black curves which pass through the focus at angles inclined ±45° from 
the horizontal direction indicate the propagation of acoustic rays which, when reaching the 
surface (horizontal line) define the inner and outer radii of the pupil. The colored curves 
indicate wavefronts (separated by a wavelength) of sample components of the egression (red) 
and ingression (blue) for this focus depth. The wavelength of the wave components shown 
here is 5.9 Mm, corresponding to p-modes with a temporal frequency of 5 mHz. If set in 
motion, the red wavefronts would diverge away from the focus, while the blue wavefronts 
would converge towards the focus. 
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Fig. 4. — Maps of the p-mode travel-time perturbations determined from lateral-vantage 
HH using the east/west quadrant pair and a focus depth of 0.7 Mm below the simulation 
surface for the following frequency bandpasses: a) 2.5 - 3.5 mHz, b) 3.5-4.5 mHz, c) 4.5-5.5 
mHz, d) 5.5-6.5 mHz, e) the average of all four frequency bandpasses, f) the time average of 
the x component of the horizontal velocity, v x , at a depth of 0.7 Mm. The color bar gives the 
units of the travel-time perturbation for panels a-e, while the values for the velocity (panel 
f) range from 1.1 km/s (black) to -1.1 km/s (white). A positive value of v x corresponds to 
a flow towards the right. 
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Fig. 5. — Comparisons of the observed p-mode travel-time perturbations determined from 
lateral- vantage HH with the perturbations computed from a forward model for several focus 
depths: a) the observed, b) model, and c) residual (observed - model) perturbations at a 
focus depth of 3 Mm, d-f) the corresponding perturbations at a focus depth of 5 Mm, g-i) 
the corresponding perturbations at a focus depth of 7 Mm below the surface 
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Fig. 6. — Measurements of the deviation from the forward models of the travel-time pertur- 
bations as functions of the degree of spatial smoothing of the observed perturbations. The 
FWHM indicates the full-width-at-half-maximum of the two dimensional Gaussian used to 
smooth the travel-time maps. The dashed lines indicate the results using 1-mHz wide fre- 
quency filters centered at 3 mHz for focus depths of 3 Mm (green), 5 Mm (blue), and 7 Mm 
(red) below the surface. The solid lines show the results for the frequency-averaged mea- 
surements. The diamonds indicate the standard deviations of the travel-time perturbations 
across the four independent frequency filters. The X symbols indicate similar standard devi- 
ations (across only three filters) for MDI observations (see text). The colors of the symbols 
indicate the same focus depths as the lines. 
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Fig. 7. — Maps of frequency-averaged travel-time perturbations measured from the simula- 
tions and from a cropped region of the Sun observed with SOHO/MD1. The top three panels 
show the results for the simulations at focus depths of a) 3 Mm, b) 5 Mm, and c) 7 Mm, and 
the bottom three panels show the same depths for the MDI observations. Only frequency 
bandpasses centered at 3, 4, and 5 mHz were used to make the simulation and MDI averages 
displayed here. 
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Fig. 8. — a) a vertical slice at y — 24 Mm through the time-averaged x component of 
the flow from the simulation, b) a slice through the travel-time kernel for a focus depth of 
3 Mm, the color shows values in units of s Mm -3 /(m/s). The travel-time kernel is defined 
as K/27rz/o, where K is given by eq. [I] for the antisymmetric phase shift computed between 
the east and west pupil quadrants, and vq = 3 mHz. The travel-time kernel is negative and 
relates a positive (negative) flow in the x direction to a travel-time decrease (increase), c) 
a vertical slice, at y — 24 Mm, through the horizontal convolution of the kernel with the 
x component of the flow field; the color bar has units of s Mm -1 , d) the total travel-time 
perturbation (blue), the target signal (red) and the remaining (total minus target) signal 
(black) as functions of the x coordinate of the focus position. The target signal is integrated 
over a range of depth from 2 to 6 Mm. 
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Fig. 9. — Same as Figure M except now the focus depth is 5 Mm. The target signal is 
integrated over a range of depth from 4 to 8 Mm. 
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Fig. 10.— Same as figure Figure [8] except now the focus depth is 7 Mm. The target 
signal is integrated over a range of depth from 5 to 10 Mm. In panel c, at x ~ 30 Mm 
notice the positive contribution to the travel-time perturbation from the near-surface. This 
contribution is due to flows located in the near-surface lobes of the kernel and with opposite 
sign to the flow at the focus point. 



